Mott transition and suppression of orbital fluctuations in orthorhombic 3d 1 perovskites 
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Using t2 g Wannier-functions, a low-energy Hamiltonian is derived for orthorhombic 3d 1 transition- 
metal oxides. Electronic correlations are treated with a new implementation of dynamical mean-field 
theory for non-cubic systems. Good agreement with photoemission data is obtained. The interplay 
of correlation effects and cation covalency (GdFeC>3-type distortions) is found to suppress orbital 
fluctuations in LaTiOa , and even more in YTiOs , and to favor the transition to the insulating state. 
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Transition-metal perovskites have attracted much in- 
terest because of their unusual electronic and magnetic 
properties arising from narrow 3c? bands and strong 
Coulomb correlations . The 3d 1 perovskites are partic- 
ularly interesting, since seemingly similar materials have 
very different electronic properties: SrV0 3 and CaV0 3 
are correlated metals with mass-enhancements of respec- 
tively 2.7 and 3.6 while LaTi0 3 and YTi0 3 are Mott 
insulators with gaps of respectively 0.2 and 1 eV 0. 

In the Mott-Hubbard picture the metal-insulator tran- 
sition occurs when the ratio of the on-site Coulomb re- 
pulsion and the one-electron bandwidth exceeds a criti- 
cal value U c /W, which increases with orbital degeneracy 
0, 0] . In the AB0 3 perovskites the transition-metal ions 
(B) are on a nearly cubic (orthorhombic) lattice and 
at the centers of corner-sharing 0$ octahedra. The 3d 
band splits into pd7r-coupled t2 g bands and pdcr-coupled 
e g bands, of which the former lie lower, have less O char- 
acter, and couple less to the octahedra than the latter. 
Simplest theories for the d 1 perovskites Q are therefore 
based on a Hubbard model with 3 degenerate, ^ -filled 
t2 g bands per B-ion, and the variation of the electronic 
properties along the series is ascribed to a progressive 
reduction of W due to the increased bending of the pdir 
hopping paths (BOB bonds). 

This may not be the full explanation of the Mott tran- 
sition however, because a splitting of the ti g levels can 
effectively lower the degeneracy. In the correlated metal, 
the relevant energy scale is the reduced bandwith asso- 
ciated with quasiparticle excitations. Close to the tran- 
sition, this scale is of order ~ ZW, with Z ~ 1 — U/U c , 
and hence much smaller than the original bandwith W. 
A level splitting by merely ZW is sufficient to lower 
the effective degeneracy all the way from three-fold to a 
non-degenerate single bandQ. This makes the insulating 
state more favorable by reducing t/ c /W|a,l3- Unlike in 
e g -band perovskites, such as LaMn0 3 , where large (10%) 
cooperative Jahn- Teller (JT) distortions of the octahedra 
indicate that the orbitals are spatially ordered, in the t2 g - 
band perovskites the octahedra are almost perfect. The 
t2 g orbitals have therefore often been assumed to be de- 



generate. If that is true, it is conceivable that quantum 
fluctuations lead to an orbital liquid Q rather than or- 
bital ordering. An important experimental constraint on 
the nature of the orbital physics is the observation of an 
isotropic, small-gap spin-wave spectrum in both insula- 
tors |8j. This is remarkable because LaTi0 3 is a G-type 
antiferromagnet with 7V=140K, m—0A5fiB, and a 3% 
JT stretching along a [9j, while YTi0 3 is a ferromagnet 
with Tc=30K, m ^0.8^b, and a 3% stretching along y 
on sites 1 and 3, and x on 2 and 4 [ljj (see Fig. 1). 

We shall find that the t2 g degeneracy is lifted at the 
classical level. This is not due to the small JT distortions 
via OB pd7r-coupling, but to the GdFe0 3 -type distortion 
which tilts the corner-sharing octahedra around the b- 
axis (by 0,9,12, and 20°) and rotates them around the 
c-axis (by 0,7,9, and 13°), as we progress from cubic 
SrV0 3 via CaV0 3 and LaTi0 3 to YTi0 3 0, EH E! G3 . 
This distortion is driven by the increasing oxygen-cation 
(OA) pder-covalency, and it primarily pulls closer 4 of the 
12 oxygens neighboring a given cation^H- Moreover, 2 
to 4 of the 8 cations neighboring a given B ion are pulled 
closer Q. The t2 g orbitals couple to the OA distor- 
tion via oxygen (BOA dpn-pda), and they couple directly 
(AB dda) to the AB distortion. As seen in Fig. 1, the 
orthorhombic GdFe0 3 -type distortion also leads to qua- 
drupling of the cell. These findings are consistent with 
conclusions drawn in the most recent model Hartree-Fock 
study |l5|. The correct magnetic orders in LaTi0 3 and 
YTi0 3 were also obtained with the LDA+ U method [lfij ]. 
The predicted moment and orbital order in YTi0 3 were 
confirmed by NMR 17] and neutron scattering ^3, but 
not in LaTi0 3 . These static mean- field methods are not 
appropriate for the metallic systems, however. 

In this letter, we shall (i) present a new implementa- 
tion of a many-body method [lfl Eflj ] , which allows for 
a quantitative, material-specific description of both the 
Mott transition and the orbital physics, and (ii) use it to 
explain why some of the d 1 perovskites are metallic and 
others are insulators, why the metals have different mass 
enhancements and the insulators different gaps. Such 
properties can be described by a low-energy, multi-band 
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FIG. 1: Pbnm primitive cell (right), subcell 1 (left), and 
the occupied t^g orbitals for LaTiOa (top) and YTiOs 
(bottom) according to the LDA+DMFT calculation. The 
oxygens are violet, the octahedra yellow and the cations 
orange. In the global, cubic zyz-system directed ap- 
proximately along the Ti-O bonds, the orthorhombic 
translations are a= (1, — 1, 0) (1 + a) , b= (1, 1, 0) (1 + 0) , 
and c= (0,0, 2) (1 + 7) , with a, f3, 7 small. The Ti sites 1 
to 4 are: a/2, b/2, (a+c)/2, and (b+c)/2. The La(Y) 
afo-plane is a mirror (z <-> —2), and so is the Ti fee-plane 
(x <-> y) when combined with the translation (b-a)/2. See: 
http://www.mpi-stuttgart.mpg.de/andersen/cm/0309102.html 



Hubbard Hamiltonian 
H = H LDA 
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where n ima = a^ mo a ima , and af ma creates an elec- 
tron with spin a in a localized orbital m at site i. 
This Hamiltonian depends on how the imcr-orbitals 
are chosen. fj LDA [ s the one-electron part given by 
density- functional theory (LDA), which should provide 
the proper material dependence. Recently it has become 
feasible to solve Q u sing the dynamical mean-field ap- 
proximati on ( DMFT) [19( and to obtain realistic physical 
properties 20]. In the original LDA+DMFT implemen- 
tations it was assumed that the on-site block(s) of the 
single-particle Green function is diagonal in the space of 
the correlated orbitals, and these were taken as orthonor- 
mal LMTOs approximated by truncated and renormal- 
ized partial waves. Although these approximations are 
good for cubic t 2g systems such as SrVC>3 |2l| , they dete- 
riorate with the degree of distortion. Our new implemen- 
tation of LDA+DMFT uses a set of localized Wannicr 
functions in order to construct a realistic Hamiltonian 
Q), which is then solved by DMFT, including the non- 
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FIG. 2: t2g LDA DOS matrix (states/eV/spin/band) in the 
Wannier representation. On-site-1 elements: N XZtXZ (red) 



"yz, yz (green) 
Nxz,xy (green), 



and N xy , xy (blue). Insets: 
and N X y, yz (blue), ef = 0. 



(red), 



diagonal part of the on-site self-energy. 

For an isolated set of bands, a set of Wannier functions 
constitutes a complete, orthonormal set of orbitals with 
one orbital per band. For the d 1 perovskites we take 
the correlated orbitals to be three localized ti g Wannier- 
orbitals, and in H LDA we neglect the degrees of free- 
dom from all other bands. In order to be complete, such 
a Wannier orbital must have a tail with e.g. Opir and 
A d characters. Our Wannicr orbitals are symmetrically 
orthonormalized Nth-order muffin-tin orbitals (NMTOs) 
|22| . which have all partial waves other than Bxy, yz, 
andzx downfolded. Such a t 2g NMTO can have on-site 
e g character, and that allows the orbital to orient itself 
after the surroundings, although xy,yz, and zx refers to 
the global cubic axes defined in Fig. 1. Fourier transfor- 
mation of the orthonormalized 12x12 NMTO Hamilto- 
nian, H LDA (k) , yields on-site blocks and hopping inte- 
grals. For the on-site Coulomb interactions in Eq. (1), 
we use the common assumption that, as in the isotropic 
case, they can be expressed in terms of two parameters: 
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Ref.|2J, J=0.68eV for the vanadates and 0.64 eV for the 
titanates. Since our Hamiltonian involves only correlated 
orbitals, so that the number of correlated electrons is 
fixed, the double-counting correction amounts to an ir- 
relevant shift of the chemical potential. H is now solved 
within DMFT, i.e. under the assumption that the com- 
ponents of the self-energy between different sites can be 
neglected. As a result, the self-energy can be obtained 
from the solution of an effective local impurity model 
which involves only 3 correlated orbitals. In contrast 
to all previous studies, we take all components of the 
self-energy matrix S mm ' between different Wannier func- 
tions on a given B-site into account |2^| . From this 3x3 
matrix, by use of the Pbnm symmetry (Fig. 1), we con- 



struct a 12x12 block-diagonal self-energy matrix. The 
latter is then used together with H LDA (k) to obtain the 
Green function at a given k-point. Fourier transforma- 
tion over the entire Brillouin zone yields the local Green 
function associated with a primitive cell and its 3x3 on- 
site block is used in the DMFT self-consistency condition 
in the usual manner. The 3-orbital impurity problem 
is solved by a numerically exact quantum Monte Carlo 
scheme [i^. To access temperatures down to 770 K, we 
use up to 100 slices in imaginary time. 10 6 QMC sweeps 
and 15-20 DMFT iterations suffice to reach convergence. 
Finally, the spectral function is obtained using the max- 
imum entropy method . 

We now present the LDA results for the four per- 
ovskites. Fig. 2 displays the on-site DOS matrix N mm > (e) 
in the representation of the xy, yz, and zx Wannier func- 
tions. SrV03 is cubic and its ti g band with a width 
I / F=2.8 eV consists of 3 non-interacting subbands, each of 
which is nearly 2D and gives rise to a nearly logarithmic 
DOS peak. In CaV03, W is reduced to 2.4 eV because 
the Wannier orbitals are misaligned by the GdFe03-type 
distortion and because some of their O 2p character is 
stolen by the increased OA covalency, which drives this 
distortion. The energy of the xy Wannier orbital (the 
center of gravity of N XVtXy ) is 80meV lower than that of 
the degenerate xz and yz orbitals, and small off-diagonal 
DOS elements appear. Going from the vanadates to 
the titanates, the effects of OA and (Ad)(B*2g) cova- 
lency increase dramatically, because now A and B are 
1st rather than 3rd-nearest neighbors in the periodic ta- 
ble. As consequences, the increased misalignment and 
loss of oxygen character reduces the bandwidths to 2.1 
and 2.0 eV in LaTiOa and YTi03, and weak hybridiza- 
tion with the Ad bands deforms the ti g band. A pseudo- 
gap which starts out as a splitting of the van Hove peak 
in CaVOs, deepens and moves to lower occupancy as 
we progress to LaTiOs and YTi03. The deep pseudo- 
gap in YTiOs is mainly caused by the hybridization with 
the Y xy orbital. The xy, yz, and zx Wannier orbitals are 
now strongly mixed and diagonalization of the on-site 
blocks of H LDA yields three singly-degenerate levels with 
the middle (highest) being 140 (200) meV above the 
lowest in LaTi0 3 , and 200 (330) meV in YTi0 3 . This 
splitting is not only an order of magnitude smaller than 
the t2 g bandwidth, but also smaller than the subband- 
widths, in particular for LaTiOs. As a consequence, the 
eigenfunction for the lowest level is occupied by merely 
0.45 electron in LaTiOs and 0.50 in YTiOs, while the re- 
maining 0.55 (0.50) electron occupies the two other eigen- 
functions. The eigenfunction on site 1 with the lowest en- 
ergy is 0.604|a;y) + 0.353|xz) + 0.7U\yz) in LaTi0 3 and 
0.619|a;y)-0.073|xz)+0.782|yz) in YTi0 3 . The splittings 
are large compared with the spin-orbit splitting (20mcV) 
and kT, and they are not caused by the JT distortions, as 
we have verified by turning them off in the calculations. 

Next, we turn to the LDA+DMFT results. Calcula- 
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FIG. 3: DMFT spectral function at T = 77QK (thick line) 
and LDA DOS (thin line). fi = 0. 



tions were performed for several values of U between 3 
and 6 eV. We found that the critical ratio U c /W decreases 
when going along the series: SrV03, CaVOs, LaTiOs, 
and YTiOs. This is consistent with the increasing split- 
ting of the t2 g levels and indicates that the Mott transi- 
tion in the d 1 series is driven as much by the decrease 
of effective degeneracy as by the reduction of bandwidth. 
The main features of the photocmission spectra for all 
four materials, as well as the correct values of the Mott- 
Hubbard gap for the insulators are reproduced by 
taking U constant ^5eV. This is satisfying, because U 
is expected to be similar for vanadates and titanates, 
although slightly smaller for the latter (2^ . In Fig. 3 
we show the DMFT spectral functions together with the 
LDA total DOS. For cubic SrVOs we reproduce the re- 
sults of previous calculations [2lll28| : the lower Hubbard 
band (LHB) is around — 1.8 eV and the upper Hubbard 
band (UHB) around 3eV. Going to CaV03, the quasi- 
particle peak looses weight to the LHB, which remains at 
— 1.8eV, while the UHB moves down to 2.5 eV. These re- 
sults are in good agreement with photoemission data [29j 
and show that SrVOs and CaV03 are rather similar, with 
the latter slightly more correlated. Similar conclusions 
were drawn in Ref.l2ll From the linear regime of the 
self-energy at small Matsubara frequencies we estimate 
the quasi-particle weight to be Z=0A5 for SrVOs and 
0.29 for CaVOs. For a k-independent self-energy, as as- 
sumed in DMFT, this yields m*/m = l/Z = 2.2 for 
SrVOs and 3.5 for CaV03, in reasonable agreement with 
the optical-conductivity values 2.7 and 3.6 |2j. 

For LaTi0 3 and YTi0 3 the LHB is around -1.5 eV, 
in accord with photoemission [3(J ■ Despite very sim- 
ilar bandwidths, the gaps are very different, 0.3 and 
1 eV, and this also agrees with experiments |3j . This 
is consistent with our findings that the ^g-level split- 
tings are smaller and (U — 2 J) c /W is larger in LaTiOs 
than in YTi03, where the orbital degeneracy is effec- 
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tively 1. Diagonalization of the matrix of occupation 
numbers reveals that for the titanates one orbital per site 
is nearly full, in contrast to LDA. It contains 0.88 electron 
in LaTi03 and 0.96 in YTi03. The orbital polarization 
increases around the metal-insulator transition and be- 
comes complete thereafter. Thus, for the vanadates, each 
orbital is approximately 1/3 occupied for all U in the 
range Oto 6 eV. The nearly complete orbital polarization 
found by LDA+DMFT for the two insulators indicates 
that correlation effects in the paramagnetic Mott insu- 
lating state considerably decrease orbital fluctuations, 
and makes it unlikely that YTiC>3 is a realization of an 
orbital liquid Q • In LaTiC>3 some orbital fluctuations 
are still active, although quite weak. The occupied or- 
bital in LDA+DMFT is 0.586\xy) +0.27 5\xz) + 0.7 62\yz) 
for LaTi0 3 and 0.622|a;y) - 0.029|xz) + 0.782|yz) for 
YTiC>3. Hence, it is nearly identical with the ones we 
obtained from the LDA as having the lowest energy. 
For YTiC>3 our orbital is similar to the one obtained 
in Ref. |lj| and for LaTiC>3 it is similar to the one ob- 
tained in Refs. 0, 0|. Our accurate Wannier func- 
tions show why these orbitals (Fig. 1, left) have the 
lowest energy: The positive (blue) lobes have bonding 
3zf n — 1 = (\xy) + \xz) + \yz)) /y/3 character on the 
nearest cations -those along [111]- and the negative (red) 
lobes have bonding xy character on the next-nearest 
cations -those along [1-11]- whose oxygen surrounding- 
is favorable for this type of bond, i.e. where (Op)(Y xy) 
pdcr-hybridization is strong. The former type of AB co- 
valency dominates in LaTiC>3, while the latter dominates 
in YTiC>3, where the shortest YO bond is merely 10% 
longer than the TiO bond. The difference seen (Fig. 1, 
right side) between the orbital orders in the two com- 
pounds is therefore quantitative rather than qualitative; 
it merely reflects the extent to which the orbital has 
the 6c-plane as mirror. The two different JT distortions 
of the oxygen square is a reaction to, rather than the 
cause of the difference in the orbital orders. This differ- 
ence is reflected in the hopping integrals between nearest 
neighbors: ^=^=99 (38) meV and £ z = 105 (48) meV 
for LaTiC>3 (YTiC^). These hoppings are fairly isotropic 
and twice larger in LaTiC>3 than in YTi03. Moreover, 
the hoppings to the two excited orbitals are stronger in 
YTiC>3 than in LaTiC>3. All of this is consistent with 
LaTiC>3 being G-type anti- and YTiC>3 ferromagnetic at 
low temperature, and it warrants detailed future calcu- 
lations of the spin-wave spectra. 

In conclusion, we have extended the LDA+DMFT ap- 
proach to the non-cubic case using ab-initio downfolding 
in order to obtain a low-energy Wannier Hamiltonian. 
Applying this method to the Mott transition in 3d 1 per- 
ovskites, we have explained the photoemission spectra 
and the values of the Mott gap without adjustable pa- 
rameters, except a single value of U. The Mott transition 
is driven by correlation effects and GdFeOvtype distor- 
tion through reduction, not only of bandwidth, but also 



of effective orbital degeneracy. Correlation effects and 
cation covalency suppress orbital fluctuations in the high- 
temperature paramagnetic insulating phase of LaTiC>3 
and YTi0 3 . 

We thank J. Nuss, G. Khaliullin, E. Koch, J. Merino, 
M. Rozenberg, A. Bringer, M. Imada for useful discussions 
and the KITP Santa Barbara for hospitality and sup- 
port (NSF Grant PHY99-07949). Calculations were per- 
formed at MPI-FKF Stuttgart and IDRIS Orsay (project 
No. 021393). S.B. acknowledges support from the CNRS 
and the EU (Contract No. HPMF-CT-2000-00658). 



[1] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998).' 
[2] H. Makino et al, Phys. Rev. B 58, 4384 (1998). 
[3] Y. Okimoto, et al, Phys. Rev. B 51, 9581 (1995). 
[4] M.J. Rozenberg, Phys. Rev. B 55, R4855 (1997). 
[5] E. Koch, O.Gunnarsson, R.M. Martin, Phys. Rev. B 60, 
15714 (1999); S. Florens et al., ibid. 66, 205102 (2002). 
[6] N. Manirii, G.E. Santoro, A. Dal Corso, and E. Tosatti, 

Phys. Rev. B 66, 115107 (2002). 
[7] G. Khaliullin, S.Maekawa, Phys. Rev. Lett. 85, 3950 (2000). 
[8] B. Keimer et al, Phys. Rev. Lett. 85, 3946 (2000); 

C. Ulrich et al, Phys. Rev Lett. 89, 167202 (2002). 
[9] M. Cwik et al, cond-mat/0302087 
[10] M. Eitel et al, J. Less-Uommon Met. 116, 95 (1986). 
[11] M.J. Rey et al, J. Solid State Chem. 86, 101 (1990). 
[12] M.H. Jung, and H. Nakotte, unpublished. 
[13] By 0, 10, 13, and 21% of the average OA distance. 
[14] By 0, 3, 4, and 10% of the average AB distance. 
[15] M.Mochizuki, and M. Imada, Phys. Rev. Lett. 91, 167203 
(2003); T. Mizokawa, D.I. Khomskii, and GA. Sawatzky, 
Phys. Rev. B 60, 7309 (1999). 
[16] I. Solovyev, N. Hamada, and K. Terakura, Phys. Rev B 
53, 7158 (1996); H. Sawada, and K. Terakura, ibid. 58, 
6831 (1998). 

[17] M.Itho et al. J. Phys. Soc. Jap. 68, 2783 (1999). 

[18] J. Akimitsu et al. J. Phys. Soc. Japan 70, 3475 (2001). 

[l9] A. Georges, G. Kotliar, W. Kraut, M.J. Rozenberg, Rev. 

Mod. Phys. 68, 13 (1996). 
[20] V.Anisimov et al. J J.Phys: Cond. Matt. 9, 7359 (1997); 

A. I. Lichtenstein and M. I. Katsnelson, Phys. Rev. B 57, 

6884 (1998). 
[21] LA. Nekrasov et al, cond-mat/0211508 
[22] O.K. Andersen and T. Saha-Dasgupta, Phys. Rev B 62, 

16219 (2000); Bull. Mater. Sci. 26, 19 (2003). 
[23] R. Fresard and G. Kotliar, Phys. Rev. B 56,12909 (1997). 
[24] T. Mizokawa and A. Fujimori, Phys. Rev. B 54, 5368 (1996). 
[25] E mm / (^ m ) (lo) 7^ in non-cubic systems, also in the eigen- 

representation of the density matrix. In Fig. 2, a large 

Y mm ' ( ^ m) points to a large E mm / ( ^ m) . 
[26] J.E. Hirsch and R.M. Fye, Phys. Rev. Lett. 56, 2521 (1986). 
[27] J.E. Gubernatis, M. Jarrell, R. N. Silver and D. S. Sivia, 

Phys. Rev. B 44, 6011 (1991). 
[28] A. Liebsch, Phys. Rev. Lett. 90, 096401 (2003). 
[29] K. Maiti et al, Europhys. Lett. 55, 246 (2001); A. 

Sekiyama et al, cond- mat /0206471 
[30] A. Fujimori et al, Phys. Rev. B 46, 9841 (1992); K. 

Morikawa et al, ibid. 54, 8446 (1996). 



